Synthetic magnetic fields for ultracold neutral atoms 



O 
<N 



C/3 



O 

o 



(N 

O 

o 
o 



'i 



Y.-J. Lin^, R. L. Compton^, K. Jimenez-Garcia^'^, J. V. Porto^, and I. B. Spielman^ 

^ Joint Quantum Institute, National Institute of Standards and Technology, 

and University of Maryland, Gaithersburg, Maryland, 20899, USA and 

^ Departamento de Fisica, Centro de Investigacion y Estudios Avanzados 

del Instituto Politecnico Nacional, Mexico D.F., 07360, Mexico 

(Dated: July 5, 2010) 

Neutral atomic Bose condensates and degenerate Fermi gases have been used to realize impor- 
tant many-body phenomena in their most simple and essential forms pjt3j, without many of the 
complexities usually associated with material systems. However, the charge neutrality of these 
systems presents an apparent limitation - a wide range of intriguing phenomena arise from the 
Lorentz force for charged particles in a magnetic field, such as the fractional quantum Hall states 
in two-dimensional electron systems ^I E] • The limitation can be circumvented by exploiting the 
equivalence of the Lorentz force and the Coriolis force to create synthetic magnetic fields in rotat- 
ing neutral systems. This was demonstrated by the appearance of quantized vortices in pioneering 
experiments [SHS] on rotating quantum gases, a hallmark of superfluids or superconductors in a 
magnetic field. However, because of technical issues limiting the maximum rotation velocity, the 
met ast able nature of the rotating state and the difficulty of applying stable rotating optical lattices, 
rotational approaches are not able to reach the large fields required for quantum Hall physics [10]- 
I12j . Here, we experimentally realize an optically synthesized magnetic field for ultracold neutral 
atoms, made evident from the appearance of vortices in our Bose-Einstein condensate. Our approach 
uses a spatially-dependent optical coupling between internal states of the atoms, yielding a Berry's 
phase [13 sufficient to create large synthetic magnetic fields, and is not subject to the limitations of 
rotating systems; with a suitable lattice configuration, it should be possible to reach the quantum 
Hall regime, potentially enabling studies of topological quantum computation. 



In classical electromagnetism, the Lorentz force for a 
particle of charge q moving with velocity v in a magnetic 
field B is V X gB. In the Hamiltonian formulation of 
quantum mechanics, where potentials play a more cen- 
tral role than fields, the single-particle Hamiltonian is 
1-L = h^i)^ — qA/h)'^ /2m^ where A is the vector potential 
giving rise to the field B = V x A, fik is the canoni- 
cal momentum and m is the mass. In both formalisms, 
only the products gB and qA are important. To generate 
a synthetic magnetic field B* for neutral atoms, we en- 
gineer a Hamiltonian with a spatially dependent vector 
potential A* producing B* = V x A*. 

The quantum mechanical phase is the relevant and sig- 
nificant quantity for charged particles in magnetic fields. 
A particle of charge q traveling along a closed loop ac- 
quires a phase (j) = 27r$B/^o due to the presence of mag- 
netic field B, where ^b is the enclosed magnetic flux and 
$0 = h/q is the flux quantum. A similar path-dependent 
phase, the Berry's phase[T3], is the geometric phase ac- 
quired by a slowly moving particle adiabatically travers- 
ing a closed path in a Hamiltonian with position depen- 
dent parameters. The Berry's phase depends only on the 
geometry of the parameters along the path, and is dis- 
tinct from the dynamic contribution to the phase which 
depends upon the speed of the motion. 

The close analogy with the Berry's phase implies that 
properly designed position-dependent Hamiltonians for 
neutral particles can simulate the effect of magnetic fields 
on charged particles. We create such a spatially- varying 
Hamiltonian for ultracold atoms by dressing them in 
an optical field that couples different spin states. The 
appropriate spatial dependence can originate from the 



laser beams' profile [TO^, TT, [15] or, as here, a spatially- 
dependent laser- atom detuning [16]. An advantage of this 
optical approach compared to rotating gases is that the 
synthetic field exists at rest in the lab frame, allowing all 
trapping potentials to be time-independent. 

The large synthetic magnetic fields accessible by this 
approach make possible the study of unexplored bosonic 
quantum-Hall states, labeled by the filling factor v = 
^b/^o^ the ratio of atom number to the number of 
flux quanta. The most outstanding open questions in 
quantum-Hall physics center on states whose elemen- 
tary quasiparticle excitations are anyons: neither bosons 
nor fermions. In some cases these anyons may be non- 
ahelian^ meaning that moving them about each other 
can implement quantum gates, thus non-abelian anyons 
are of great interest for this "topological" quantum 
computation [17]. In electronic systems, the observed 

V = 5/2 quantum-Hall state may be such a system, but 
its true nature is still uncertain [18]. In contrast, the 

V = 1 bosonic quantum-Hall state with contact inter- 
actions has the same non-abelian anyonic excitations as 
the V — 5/2 state in electronic systems is hoped to [19]. 

To engineer a vector potential A* = A*f , we illumi- 
nate a ^^Rb BEC with a pair of Raman laser beams with 
momentum difference along x (Fig. la). These couple 
together the three spin states, m^ = and ±1, of the 
55'i/2,F = 1 electronic ground state (Fig. lb), produc- 
ing three dressed states whose energy-momentum disper- 
sion relations Ej{kx) are experimentally tunable. Ex- 
ample dispersions are illustrated in Fig. Ic. The lowest 
of these, with minimum at /Cmin, corresponds to a term 
in the Hamiltonian associated with the motion along x. 
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FIG. 1: Experiment summary for synthesizing magnetic fields, a, The BEG is in a crossed dipole trap in a magnetic 
field B = (Bo — b'y)y. Two Raman beams propagating along y ^ x (linearly polarized along y zL x) have frequencies cjl, 
ujl + Acjl. b, Raman coupling scheme within the F = 1 manifold: uz and e are the linear and quadratic Zeeman shifts, 
and S is the Raman detuning, c, Energy-momentum dispersion relations. The grey curves represent the states absent Raman 
coupling; the three colored curves are Ej{kx) of the dressed states. The arrow indicates the minimum at /cmin- d, Vector 
potential q" A% = /i/cmin versus Raman detuning 8. The insets show the dispersion Ei{kx) for M = (top inset) and —2El 
(bottom inset). e,f, Dressed BEG imaged after a 25.1 ms TOP without (e) and with (f) a gradient. The spin components, 
TTiF = 0, ±1, separate along y due to the Stern- Gerlach effect. 



H* ^ h\k, - fc„,in)V2m* = h\k, - g*A*/n)V2m*, 
where A* is an engineered vector potential that depends 
on an externally controlled Zeeman shift for the atom 
with a synthetic charge g*, and m* is the effective mass 
along X. To produce the desired spatially-dependent 
Al{y) (Fig. Id), generating -5*i = V x A*, we apply 
a Zeeman shift that varies linearly along y. The result- 
ing 5* is approximately uniform near y = 0, at which 
point A* = B^'y. (Here, the microscopic origin of the 
synthetic Lorentz force [20 is optical along x, depending 
upon the velocity along y; the force along y is magnetic, 
depending upon the x velocity) By these means, we en- 
gineer a Hamiltonian for ultracold atoms, that explicitly 
contains a synthetic magnetic field, with vortices in the 
ground state of a BEG. This is distinctly different from 
all existing experiments, where vortices are generated by 
phase imprinting [21, 22 , rotation[7-9 , or a combination 
thereof [23] . Each of these earlier works presents a differ- 
ent means to impart angular momentum to the system 
yielding rotation. Fig. le shows an experimental image of 



the atoms with 5* = 0. Fig. If, with B* > 0, shows vor- 
tices. This demonstrates an observation of an optically 
induced synthetic magnetic field. 

We create a ^^Rb BEG in a 1064 nm crossed dipole 
trap, loaded into the lowest energy dressed state [24 with 
atom number N up to 2.5 x 10^, and a Zeeman shift 
ujzf^Tr = gjiBB/h ~ 2.71 MHz, produced by a real 
magnetic bias field By. The A = 801.7 nm Raman 
beams propagate along y ^ x and differ in frequency 
by a constant Acjl — ^z, where a small Raman de- 
tuning 5 = /S^uol — (jOz largely determines the vector 
potential A*. The scalar light shift from the Raman 
beams, combined with the dipole trap gives an approxi- 
mately symmetric three-dimensional potential, with fre- 
quencies fxifyifz ^ 70 Hz. Here, hkL = h/{-\/2\) and 
El = h^k\/2m are the appropriate units for the momen- 
tum and energy. 

The spin and momentum states Imp^kx) coupled by 
the Raman beams can be grouped into families of states 
labeled by the momentum hkx. Each family "^{kx) = 
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FIG. 2: Appearance of vortices at different detuning gradients. Data was taken for N = 1.4 x 10^ atoms at hold time 
th = 0.57 s. a-f, Images of the \mF = 0) component of the dressed state after a 25.1 ms TOF with detuning gradient 6 from 
to 0.43 kHz /im~^ at Raman couphng HQr = 8.20 El- g, Vortex number Ny versus 6 at HQr = 5.85 El (blue circles), and 
8.20 El (red circles). Each data point is averaged over > 20 experimental realizations, and the uncertainties represent standard 
deviations (a). The inset displays Ny versus the synthetic magnetic flux ^b*/^o = Aq* {B*)/h in the BEG. The dashed lines 
indicate S below which vortices become energetically unfavorable according to our GPE computation, and the shaded regions 
show the 1-cr uncertainty from experimental parameters. 



{| — l^kx + 2/ci;,),|0, kx)^\ -\- l^kx — 2/ci,)} is composed of 
states that differ in linear momentum along x by ±2fi/cL, 
and are Raman-coupled with strength HQr. For each kx, 
the three dressed states are the eigenstates in the pres- 
ence of the Raman coupling, with energies Ej{kx)^^^. 
The resulting vector potential is tunable within the range 
-2kL < q*A%/h < 2kL. In addition, Ej(kx) includes a 
scalar potential V [16^. ^4*, y , and m* are functions of 
Raman coupling Q.r and detuning 5^ and for our typical 
parameters tti* ~ 2.5771, reducing f^ from ^ 70 Hz to 
^ 40 Hz. The EEC's chemical potential (i/h ^ 1 kHz 
is much smaller than the ~ /i x 10 kHz energy sepa- 
ration between dressed states, therefore the BEG only 
occupies the lowest energy dressed state. Further, it jus- 
tifies the harmonic expansion around (7*^4* /fi, valid at 
low energy. Hence, the complete single-atom Hamilto- 
nian is ?/ = H* + h^(kl + kl)/2m + V(r), where V{y) is 

the external potential including V (Qr^S). 

The dressed BEG starts in a uniform bias field B = 
Boy, at Raman resonance {S = 0), corresponding to 
A^ = 0[24]. To create a synthetic field 5*, we apply 
a field gradient b such that B = {Bq — b y)y^ ramp- 
ing in 0.3 s from b = to a variable value up to 
0.055 Tm~^, and then hold for th to allow the system 
to equilibrate. The detuning gradient 5 = g/J^sb /h gen- 
erates a spatial gradient in A*. For the detuning range 
in our experiment, dA'^/dS is approximately constant, 
leading to an approximately uniform synthetic field 5* 
given by 5* = dA^/dy = 6 dAl/d5 (Fig. Id). To probe 
the dressed state, we switch off the dipole trap and the 
Raman beams in less than 1 /is, projecting each atom 
into spin and momentum components. We absorption- 



image the atoms after a time- of- flight (TOF) ranging 
from 10.1 ms to 30.1 ms (Fig. lef). 

For a dilute BEG in low synthetic fields, we expect to 
observe vortices. In this regime, the BEG is described by 
a macroscopic wave function ?/;(r) = \'il){Y)\e^^^^\ which 
obeys the Gross-Pitaevskii equation (GPE). The phase (j) 
winds by 27r around each vortex, with amplitude l?/^! = 
at the vortex center. The magnetic flux ^b* results in 
Ny vortices and for an infinite, zero temperature sys- 
tem, the vortices are arrayed in a lattice [25 with den- 
sity q*B^/h. For finite systems vortices are energeti- 
cally less favorable, and their areal density is below this 
asymptotic value, decreasing to zero at a critical field 
5*. For a cylindrically symmetric BEG, B^ is given by 
(7*5*//i = 5/(27ri?2)ln(0.67i?/O where R is the Thomas- 
Fermi radius and ^ is the healing length [26^. B* is larger 
for smaller systems. For our non-cylindrically symmetric 
system, we numerically solve the GPE to determine 5* 
for our experimental parameters (see Methods). 

For synthetic fields greater than the critical value, we 
observe vortices that enter the condensate and reach an 
equilibrium vortex number Ny after « 0.5 s. Due to a 
shear force along x when the Raman beams are turned 
off, the nearly-symmetric insitu atom cloud tilts during 
TOF. While the vortices' positions may rearrange, any 
initial order is not lost. During the time of our experi- 
ment, the vortices do not form a lattice and the positions 
of the vortices are irreproducible between different ex- 
perimental realizations, consistent with our GPE simula- 
tions. We measure Ny as a function of detuning gradient 
S at two couplings, HQr = 5.85 El and 8.20 El (Fig. 2). 
For each Qr^ vortices appear above a minimum gradient 
when the corresponding field (5*) = 5 (dA^/dS) exceeds 
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FIG. 3: Vortex formation, a-f images of the \mF = 0) component of the dressed state after a 30.1 ms TOF for hold times 
th between -0.019 s and 2.2 s. The detuning gradient S /27t is ramped to 0.31 kHz /im~^ at the coupUng HQr = 5.85 El- g, 
top: time sequence of S . bottom: vortex number Ny (sohd symbols) and atom number N (open symbols) versus th with a 
population lifetime of 1.4(2) s. The number in parenthesis is the uncorrelated combination of statistical and systematic 1-cr 
uncertainties. 



the critical field B^. (For our coupling B^ is only approx- 
imately uniform over the system and (5*) is the field av- 
eraged over the area of the BEC.) The inset shows Ny for 
both values of Qr plotted versus ^b*/^o = Aq*{B*)/h, 
the vortex number for a system of area A = TvRxRy 
with the asymptotic vortex density, where Rx (Ry) is 
the Thomas- Fermi radius along x (or y). Since the sys- 
tem size, and thus 5J, are approximately independent 
of Qr, we expect this plot to be nearly independent of 
Raman coupling. Indeed, the data for HQr = 5.S5El 
and 8.20£^i, only deviate for A^^ < 5, likely due to the 
intricate dynamics of vortex nucleation[27 . 

Figure 3 illustrates a progression of images showing 
vortices nucleate at the system's edge, fully enter to an 
equilibrium density and then decay along with the atom 
number. The time scale for vortex nucleation depends 
weakly on 5*, and is more rapid for larger 5* with more 
vortices: It is ~ 0.3 s for vortex number Ny > S^ and in- 
creases to ^ 0.5 s for Ny = 3. For Ny = 1 {B* near 5*), 
the single vortex always remains near the edge of the 
BEC. In the dressed state, spontaneous emission from 
the Raman beams removes atoms from the trap, caus- 
ing the population to decay with a 1.4(2) s lifetime, and 
the equilibrium vortex number decreases along with the 
BEC's area. 

To verify the dressed state has reached equilibrium, 
we prepare nominally identical systems in two different 
ways. First, we vary the initial atom number and mea- 
sure Ny as a function of atom number A/" at a fixed hold 
time th = 0.57 s. Second, starting with a large atom 
number, we measure both Ny and N, as they decrease 



with th (Fig. 3). Figure 4 compares Ny versus TV mea- 
sured with both methods, each at two detuning gradients 
corresponding to fields B^ < ^2 • The data show Ny as a 
function of N is the same for these preparation methods, 
providing evidence that for t^ > 0.57 s, Ny has reached 
equilibrium. As the atom number A" falls, the last vor- 
tex departs the system when the critical field - increasing 
with decreasing A' - surpasses the actual field. 

In conclusion, we have demonstrated optically synthe- 
sized magnetic fields for neutral atoms resulting from the 
Berry's phase, a fundamental concept in physics. This 
novel approach differs from experiments with rotating 
gases, where it is difficult to add optical lattices and ro- 
tation is limited by heating, metastability, and the diffi- 
culties to add large angular momentum, preventing ac- 
cess to the quantum-Hall regime. A standout feature in 
our approach is the ease to add optical lattices. For ex- 
ample, the addition of a 2D lattice makes it immediately 
feasible to study the fractal energy levels of the Hofs- 
tadter butterfly [28 . Further, a ID lattice can divide the 
BEC into an array of 2D systems normal to the field. A 
suitable lattice configuration allows access to the u ^ 1 
quantum-Hall regime, with an ensemble of 2D systems 
each with ^ 200 atoms, and with a realistic ~ /c^ x 20 nK 
interaction energy. 
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FIG. 4: Equilibrium vortex number. Vortex number 
Nv versus atom number N at detuning gradient Si/27t = 
0.26 kHz /im~^ (red circles) and S2/27T = 0.31 kHz /im~^ 
(black circles), corresponding to synthetic fields Bi < B2, at 
Raman coupling HQr = 5.85 El. The two data points with 
the largest N show representative 1-cr uncertainties, estimated 
from data in Fig. 2g. We vary N by its initial value with a 
fixed hold time th = 0.57 s (solid symbols), and by th with a 
fixed initial N (open symbols) . The vertical dashed lines indi- 
cate N below which vortices become energetically unfavorable 
computed using our GPE simulation. The shaded regions re- 
flect the l-cr uncertainties from the experimental parameters. 
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Methods Summary 



taneously, we decrease the dipole beam power by 20%, 
producing our ^ 40 Hz trap frequency along x. Addition- 
ally, the detuning gradient S y makes the scalar potential 
V anti-trapping along ^, reducing fy from 70 Hz to 50 Hz 
for our largest S . Spontaneous emission from the Raman 
beams decreases the atom number to N ^ 2.5 x 10^ for 
th = 0, with a condensate fraction of 0.85. 



A. Dressed state preparation 

We create a ^^Rb BEC in a crossed dipole trap [29], 
with N ^ 4.7 X 10^ atoms in \F = I.tuf = -1). The 
quadratic Zeeman shift is he = 0.61^l for ujz/^tt = 
gjiBB/h ^ 2.71 MHz, where g is the Lande ^f-factor. To 
maintain J = at the EEC's center as we ramp the field 
gradient b , we change qjibBq by as much as 7 El. Simul- 



B. Numerical method 

We compare our data to a finite temperature 2D 
stochastic GPE (SGPE[30]) simulation including the 
dressed state dispersion E{kx,y) that depends on y 
through the detuning gradient S' . We evolve the time- 
dependent projected GPE 
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V projects onto a set of significantly occupied modes, 
and g2D parameterizes the 2D interaction strength. The 
SGPE models interactions between the highly-occupied 
modes described by ip and sparsely occupied thermal 
modes with dissipation and an associated noise term. We 
approximately account for the finite extent along z by 
making g2D depend on the local 2D density . For low 



temperatures this 2D model correctly recovers the 3D 
Thomas-Fermi radii, and gives the expected 2D density 
profile. These quantitative details are required to cor- 
rectly compute the critical field or number for the first 
vortex to enter the system, which are directly tied to the 
2D condensate area. 



